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Abstract 

We study the conditions for soliton-like wave propagation in the Photore- 
fractive (PR) and electro-optic (i.e., Pockels) material, by using Nonlinear 
Schrodinger (NLS) equation. The complete NLS equation is solved analyt- 
ically and numerically by transforming it into the phase space. Our results 
clearly show the existence of both the dark and bright solitary solutions for the 
PR medium. Interestingly, however, we find only one bright solitary solution 
in the Pockels case and there is no evidence of any dark solitary solution. 
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I. INTRODUCTION 



In a uniform Kerr-like nonlinear medium, the transverse profile of an optical beam prop- 
agates unchanged since the diffraction is balanced by self-focusing. Such self-focused waves 
are called solitons, and are described by the Nonlinear Schrodinger (NLS) equation in space- 
time dimension - 0j . Zakharov and Shabat have solved such NLS equation for the Kerr 
medium by using Inverse Scattering Transformation (1ST) 0. Self-focused optical beams 
in other non-Kerr, nonlinear media have also been investigated for the possibility of gen- 
erating soliton like motion. For example, in the case of an optical beam propagation in a 
Photorefractive (PR) material, charge carriers are photo-excited and redistributed around 
the beam. This charge distribution produces a space charge field, which in turn modulates 
the refractive index(the electro-optic effect) ||, |J. When such a system exhibits optical 
nonlinearity, the related modifications in the refractive index induces either a self-focusing 
or self-defocusing of the incident beam. With an appropriate choice of the experimental 
parameters, the input beam profile can be made to converge asymptotically to a soliton 
state whose PR nonlinearity compensates for the diffraction, and the beam profile remains 
unchanged as it propagates. Montemezzani and Gunter have investigated the profile of the 
one-dimensional bright spatial solitons by using a power series expansion J7[. They managed 
to calculate the first four terms of the series that describe the profile of the bright solitons. 
However, a complete temporal and spatial solution is yet to emerge. 

In this work, we first attempt to solve the complete NLS equation for the Pockels medium. 
It is extremely difficult to solve this equation with 1ST, so we invoke the technique of 
transforming the NLS equation into phase space and determine its exact solution. We then 
extend the technique to solve the NLS equation for the PR medium. In fact, we find that 
the technique of phase space transformation is so general that it enables us to solve the 
NLS equation for any polynomial potential that appear in most optics and plasma physics 
problems. 
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II. SOLITARY WAVE IN DISPERSIVE POCKELS MEDIA 



Let us consider the spatial and temporal variation of the electric field E(z, t) of an optical 
pulse propagating through a linear electro-optic medium as ||: 

E(z, t) = A(z, t)exp i(k z — u t) + c.c. (1) 

Here z and t are the space and time coordinates and ko and ujq are the linear components 
of the wavevector and frequency of the electric field. The total wavevector k can be written 
as i: 

dk 1 d 2 k 

k = k + AK NL + (w- LO )(-j^)u=u + ^(^ ~ ^o) 2 (^)^ o; ( 2 ) 

the nonlinear contribution to the propagation constant is: 

AK NL = 1 \E(z,t)\= 1 \A(z,t)\ (3) 

where 7 is Pockels constant. Defining the retarded time t a.s r = t — z/v g , the NLS equation 
can be written as ||, 0: 

dA .(3 2 d 2 A , Ar ^ ~ 

where f3 2 is a measure of group velocity dispersion (GVD) and is defined in terms of wavevec- 
tor as: 

/3 2 (3 o)u=luo ~\ ( 7 \")a>=u;o ( T~] )w=wo" (5) 

duj 1 doj Vg(uj) v l g du 

Using Eqn.(3), the NLS equation for the propagating pulses in a Pockels medium is: 

dA .(3 2 d 2 A . . ~. ~ 

_ + ! |_ = !7 |^, (6) 

The effect of GVD on propagation of pulse in a linear dispersive medium can be consid- 
ered by setting 7 = in Eqn.(6), 

dA (3 2 d 2 A 
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Eqn.(7) can be solved using the Fourier method ||, which results in the familiar result 
that the long wavelength components of an optical pulse propagate faster than the short 
wavelength components when the GVD parameter, 02, is positive, and vice versa. 

Another interesting phenomenon that occurs in the nonlinear media is self-phase modu- 
lation effect (SPM) , which can be determined in the absence of dispersion (f3 2 = 0). When 
the SPM parameter 7 is negative, the short wavelength components of an optical pulse prop- 
agate faster than the long wavelength components, and vice versa ||, |§. So when fl 2 > 
and 7 < or /9 2 < and 7 > 0, the SPM effect is canceled by GVD effect, and optical 
pulse can propagate in a medium with an unchanged shape. Let us consider a solution to 
the Eqn.(6) as follows, 

A(z,t) = F(r)exp (iTz), (8) 
where T is defined as the wave propagation constant. Defining new variables as: 

(9) 



the Eqn.(6) becomes 

d 2 F(r) 



+ sgn(X)\X\F(r) + sgn(^\F 2 (r) = 0, (10) 



dr 2 

where sgn(A) and sgn(/x) takes the value of plus or minus one. 
Making a further change of variables as follows: 

UJr = T, 



F(r) = tor), (11) 



q> _ dQ(T) 



H 

[T) 
dT ' 



Eqn.(10) reduces to: 



Q" + sgn{X)Q + sgn{n)Q 2 = 0. (12) 



Let us consider this equation as an Euler-Lagrange equation of a Hamiltonian system. 
Multiplying Eqn.(12) by Q 1 and integrating with respect to time, the corresponding Hamil- 
tonian function can be derived as the following form: 



H(Q, Q') = Q' 2 + sgn{\)Q 2 + f ^n(/i)Q 3 = h, (13) 

where h is the constant of integration and can be determined from the initial conditions. 
The Hamiltonian in Eqn.(13) represents a dynamical system and its behavior is investigated 
by considering the quantities A, /i and h as parameters. 

When the parameters A < and \i > 0, Eqn.(13) can be expressed as: 

H(Q, Q') = Q' 2 + \q z - Q 2 = h. (14) 



We can find the extreme limits (extrema points) of the Hamiltonian |L0| by setting 

(15) 



( 9h_ Q 

dQ' u ' 



— = 

We obtain h = and h = — | as the extrema points. Now we can plot Q' 2 — Q 2 + |Q 3 = h 
for different values of h, which are: h = 0, — | < h < , h > and /i = — | (see Fig. (1)). 
For h — 0, the figure suggests a hyperbolic solution. Putting ft, = in Eqn. (14) we obtain 
the following differential : 

dQ ±dT (16) 



^Q 3 + Q 2 
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by integration we can obtain the following solution 



Q = lsech 2 Q, (17) 



writing in terms of the original variables, we obtain 



A(z,t) = ^sech 2 (^{t-z/v g ))eF*. (18) 

When (t — z/vg) — > ±oo, A(z,t) — > and hence this is a bright solitary solution, a profile 
of which is shown in Fig. (2). 

In the other case, A > and /j> Owe find: 

H(Q,Q') = Q' 2 + ^Q 3 + Q 2 = h. (19) 



Calculations similar to above can be done to get the extremum point to be (Q, Q') = (0, —1). 
We plot H (Q, Q') for different values of h in Fig. (3). We have hyperbolic solution for h — |. 
For this value of h, Eqn. (19) can be integrated as before as: 

Q(T) = -^tanh 2 (|) + i (20) 

Thus, 

A(z,t) = ^[-\tanh\^{t-z/v g )) + \y^ (21) 

The profile of the solution is given in Fig. (4). Since this amplitude assumes negative 
values, it is a unphysical solution and as such no dark solitons are possible. The case of 
(A < and fi < 0) or (A > and /i < 0) also lead to non-physical solution. 

III. SOLITARY WAVE IN PR MEDIA 

The nonlinear contribution to the propagation in PR media can be written as (see ap- 
pendix) : 

AA - = 7 4Tum)' (22) 

where 7 is the electronic coupling parameter, I p is the intensity of the pump beam and 
I s (z,t) is the spatial and temporal variation of the signal beam. When this form of the 
nonlinear contribution is substituted into the NLS Eqn. (4), we obtain: 

dA (3 2 d 2 A . I ~ 

^ + t Y^ = tl J-^T s A (23) 

where /3 2 is a measure of the dispersion of the group velocity. 

To achieve the separation of variables z and r, we seek a solution of the form: 

A(z,r) = JT p F(T)exp (iFz), (24) 

where T is defined as the pulse propagation constant and F{r) is a real function of r. With 
the definition: 
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we have : 

^1 + sgn(X)\X\F(r) + sgn(^\ y-^y = 0, (26) 

With a change of variables similar to Eqn. (11), Eqn. (26) reduces to: 

Q' 2 + sgn(X)Q + sgn^)-^-^ = 0, (27) 

where a 2 = H = j^j. Using Eqn. (22) and the condition for bright soliton [11|, T, soliton 
propagation constant is equal to the PR nonlinearity AK^l, i-e. &k nl = 1 which further 
gives us the condition that |r|/|7| < 1 (since I s /I p is always positive ). Similarly for dark 
solitons we have -r-£ — = ln ( 1 + Q o) _ gj nce Mh-Qo) < ^ f or ever y va ] ue Q f g Q; we have again 

|r|/H<i- 

Similar to last section, we multiply Eqn. (27) by Q' and integrate, we obtain the corre- 
sponding Hamiltonian function: 

H(Q, Q') = Q' 2 + sgn(\)Q 2 + sgn(fi) ln(a 2 + Q 2 ) = h, (28) 

where h is a constant of integration. 

We look in the phase space for solitary solutions of the above equation. Using Eqn. (15) 
we find as before the extrema points of the Hamiltonian to be: 

h = s 9 n(X)[-^ - a 2 } + sgnfr) ln(-g^) 

{ ZJ > 

h = sgn(/J,) ln(a 2 ) 
If sgn(X) = sgn(fi) then only the second value of h is valid. 

First we consider the case of \x < 0, A > 0, where the two values of h are h = 1 — a 2 and 
h = — ln(a 2 ). It can be easily seen that there is no solitary solution for the first value of h. 
For the second value of h Eqn. (28) can be expressed as: 

H(Q, Qf) = Q' 2 + Q 2 - ln(a 2 + Q 2 ) = - In a 2 (30) 



The shapes of the orbits for different values of a 2 are shown in Fig. (5). We obtain the 
following differential : 



Integrating this equation, we find the profile of the solution as plotted in Fig. (6). 
For the case \i > 0, A < 0, the solitary wave solution is located at h = a 2 — 1 and substituting 
h in the Hamiltonian of Eqn.(28), we can derive a dark solitary solution. The phase space 
orbits are given in Fig. (7) and the profile of the dark solution is given in Fig. (8). 
For all other cases (fx < 0, A < and /x > 0, A > 0), there are no solitary solutions. 



We started from the NLS equation and solved it for a time- dependent potential in a 
medium with Pockels effect. Our analytical solutions show that there is only one possible 
case where the soliton-like wave can propagate in the Pockels material. 

In the second part, we solved the NLS equation for PR medium to find the conditions 
under which the soliton-like waves can propagate. The conditions we found are: a) /i < 
0, A > that leads to the bright solitary waves and b)/i > 0, A < 0, which leads to the dark 
solitary waves. 



dQ 



= ±dT. 



(31) 




IV. CONCLUSIONS 
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V. APPENDIX 



Nonlinear index of refraction of the PR medium be derived from the two beam coupling 
(pump and signal) inside a PR material. Nonlinear polarization that is produced by the 
coupling of two beams is given by || , |TJ| : 

p nl = (tn-r + c . c .)(A s e iks - r + A„e ikp - r ) (32) 

The part of nonlinear polarization that has the spatial variation e^ik s .r) can act as a phase- 
matched source for the signal wave as: 

NL _ p ik s -r _ - ^leffEm \A p \ 2 A s 

Ps ~ 4n Ape ~~ l 4tt \A p \ 2 + \A s f [dd) 
Using the nonlinear relationship between P^ L and e NL and Eqn.(33), we can write 

Anpf L = e NL E s = e NL A s e^ r = -ie^ eff E m -^-A s e ik - r (34) 

■»p ~r -'s 

as a reseult of which the value of e NL is equal to: 

e NL = e^ eff E m -^-. (35) 

J-p ~r -'s 

Now, given the value of e NL , we can calculate the Ak^L for the PR medium by substi- 

2fc c 2 



tuting e = n 2 and k = — in the nonlinear relation A/ctvl — -^r^^ NL as follows: 



A^atx = Ij^TT, (36) 

■*p ~r -'s 

where we have defined a parameter 7 as: 

7 = ^ 3 7e//£m- (37) 

Eqn. (36) is the potential that we use in the (NLS) equation for PR media in Eqn. (23) of 
the main text. 
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FIGURES 

FIG. 1. Phase space diagram of Q /2 -Q 2 + §Q 3 = h for different values of h. The curves 1,2, 3 
and 4 correspond to h = 0, — A < h < 0(h = — A), h > 0(/i = A) and h = — A, respectively. 



FIG. 2. Normalized amplitude profile of \A(z,t)\ = f^sec/i 2 (^J% - z/v g )) for A < and 

fi > 0, which shows bright solitary solution. 

FIG. 3. Phase space diagram of Q' 2 + Q 2 + |Q 3 = h for different value of h. The curves 1,2,3 

and 4 correspond to h = g, < h < \(h = \),h< 0(h = 0) and h > |, respectively. 

FIG. 4. Normalized amplitude profile of = j^[-ftan/i 2 (^(t - z/v fl )) + ±] for A > 

and jU > in Pockels media. This profile shows that there is no solitary solution possible under 

this condition. 

FIG. 5. Phase space diagram of Q' 2 + Q 2 - ln(a 2 + Q 2 ) = - In a 2 for different value of a 2 . The 



FIG. 6. Amplitude profile of the differential , = = = = = = ±dT solved numerically 

yJ\n{a 2 +Q z )-Q A -\n(a 2 ) 

a 2 = 0.6, which shows the bright solitary solution. 

FIG. 7. Phase space diagram of Q' 2 - Q 2 + ln(a 2 + Q 2 ) = a 2 - 1 for different value of a 2 . The 
ves 1,2,3 and 4 correspond to a 2 = 0.2, 0.4, 0.6 and 0.8, respectively 

FIG. 8. Amplitude profile of the differential -— d ? 9 ^ = ±dT solved numerically for 

y/ Q 2 -ln(a 2 +Q z )a 2 -1 
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